Scale Invariant Dynamics of Surface Growth 
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We describe in detail and extend a recently introduced 
nonperturbative renormalization group (RG) method for sur- 
face growth. The scale invariant dynamics which is the key 
ingredient of the calculation is obtained as the fixed point of 
a RG transformation relating the representation of the mi- 
croscopic process at two different coarse-grained scales. We 
review the RG calculation for systems in the Kardar-Parisi- 
Zhang universality class and compute the roughness expo- 
nent for the strong coupling phase in dimensions from 1 to 
9. Discussions of the approximations involved and possible 
improvements are also presented. Moreover, very strong ev- 
idence of the absence of a finite upper critical dimension for 
KPZ growth is presented. Finally, we apply the method to 
the linear Edwards- Wilkinson dynamics where we reproduce 
the known exact results, proving the ability of the method to 
capture qualitatively different behaviors. 



I. INTRODUCTION 

The idea that scale invariance is at the origin of criti- 
cal phenomena associated with equilibrium second order 
phase transitions has proven to be very fruitful. The 
analysis of scale transformations in equilibrium statisti- 
cal systems, now known as renormalization group (RG), 
has indeed allowed for the explicit calculation of critical 
exponents and, moreover, has led to the introduction of 
new fundamental concepts such as scaling and universal- 
ity. 

The extension of the RG approach to non-equilibrium 
phenomena, where scale invariance is widely observed, 
and the identification of new universality classes, is of 
great importance from both theoretical and practical 
points of view. Technically, the RG ideas can be im- 
plemented in different ways. The most standard one 
for systems at equilibrium is to consider their station- 
ary probability distribution written in terms of contin- 
uum coarse-grained fields, and study them perturba- 
tively around their corresponding upper critical dimen- 
sion. The most systematic way to extend the previous 
methods to non-equilibrium systems, where in general 
the stationary probability distribution is not known, is 
to cast them into a continuum dynamical equation [1|, 
or equivalently into a generating functional or action [ 2 [ . 
This last one can, in principle, be treated using the same 
perturbative techniques developed to deal with equilib- 
rium systems. However, there are some cases where per- 
turbative methods around a mean field solution are not 
suitable. In these cases the e-expansion fails to give in- 



formation on the relevant physics. This turns out to be 
governed by a strong coupling, perturbatively inaccessi- 
ble fixed point. The prototypical example of this class 
of systems is the well known Kardar-Parisi- Zhang (KPZ) 
equation for surface growth |^, where the properties of 
rough surfaces have not been so far explained satisfac- 
torily in generic spatial dimension. This is a problem 
of great theoretical importance since the KPZ describes 
not only the properties of rough surfaces 0-^, but is 
also related to the Burgers equation of turbulence , to 
directed polymers in random media and to systems 
with multiplicative noise Q. In particular, one of the 
most debated issues in this context is the existence of an 
upper critical dimension, above which the system is well 
described by the nontrivial infinite dimensional limit [ p^ . 

Although the usual approach fails for KPZ, the pres- 
ence of generic scale invariance suggests that also for this 
system the basic idea of the RG approach should be ap- 
plicable in some form. 

Real space approaches have proven useful wherever 
standard perturbative techniques fail ]rT|-|l3t. This, for 
example, is the case of fractal growth, and in partic- 
ular for Diffusion-Limited- Aggregation iQ. However, 
the attempts to apply standard real space techniques 
to the KPZ problem (and to surface growth in general) 
fail because of a fundamental technical difficulty: The 
anisotropy of the scaling properties of the system. That 
is, in order to cover with blocks (in the Kadanoff sense 
p3[) a surface, isotropic blocks cannot be used: Lengths 
in different directions must scale in different ways, and 
the relative shape of blocks has to depend upon the scale 
via an exponent that is unknown. This makes conceptu- 
ally non-trivial the application of real space RG proce- 
dures to surface growth processes. 

In this paper we investigate the scale invariant prop- 
erties of generic interface growth processes through the 
introduction of a real-space method. To achieve this goal 
we introduce some new ingredients permitting us to over- 
come the aforementioned problem. In particular, we in- 
troduce the idea that the statistical properties of growing 
surfaces on large scales can be described in terms of an ef- 
fective scale invariant dynamics for renormalized blocks. 
Such dynamics is the fixed point of a RG transforma- 
tion relating the parameters of the dynamics at different 
coarse-graining levels. The study of the RG flow, of the 
fixed points and of their stability gives the universality 
classes and their associated exponents. 

As a first application of the method, we study the KPZ 
growth dynamics and obtain accurate estimates for the 
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roughness exponent (when compared with numerical re- 
sults) in spatial dimensions from d = 1 to d = 9. Fur- 
thermore, an analytical approximation allows us to ex- 
clude the existence of a finite upper critical dimension 
for KPZ dynamics and suggests that the roughness ex- 
ponent decays as l/d for large dimensions, shedding light 
on a currently much debated issue. 

In order to show the generality of the new real space 
scheme and test its accuracy, we also apply it to the well 
known linear theory, the Edwards- Wilkinson (EW) equa- 
tion. We reproduce the expected behavior in different 
dimensions, confirming the general applicability of the 
method. 

The paper is organized as follows. In section II we 
present the general RG method, the main concepts, the 
basic equations, and discuss all the approximations in- 
volved. In section III we review some results associated 
with KPZ growth and apply the new RG method to such 
problem. We present some simple analytical approxima- 
tions, explicit results for spatial dimensions up to d = 9, 
and discuss the large dimensional limit in detail. In sec- 
tion IV we report results on the analysis of the Edwards- 
Wilkinson equation. In section V a critical discussion 
of the method and of the results is reported. Partial ac- 
counts of the work presented here have already been pub- 
lished recently, with a slightly different notation [jl^^ • 

II. REAL SPACE RG FOR SURFACE GROWTH 

In order to present the RG method let us consider 
a generic surface growth model where the height is a 
single- valued function h(x,t), with x the position in a 
d-dimensional substrate and t denoting time. The pos- 
sibility of having overhangs will not be considered here, 
as they are known to be irrelevant for the asymptotic 
behavior of KPZ-like growth The generic growth 

model under consideration can be either described at the 
microscopic level by a stochastic equation or by a discrete 
dynamical rule. In the first case h and x are continuous 
variables, while in the latter they are discrete. 

The roughness of a system, when considered on a sub- 
strate of linear size L, is defined by 

X 

where 

Ht)^]TiJ2h{x,t). (2) 

X 

If we start the growth process from a flat configuration, 
for short times the roughness grows as 

W{L,t)^t'^ (3) 

until it reaches a stationary state characterized by 



W{L) - L". (4) 

The crossover between the two behaviors occurs at a 
characteristic time tg, that scales with L as L^. This 
is the time scale over which correlations decay in the sta- 
tionary state. The exponents a, P and z are to a large 
extent universal for many different growth processes, and 
are related by the trivial scaling relation /3 = a/z. 

We now introduce the real space renormalization group 
(RSRG) procedure aimed at the study the stationary 
state and in particular at the determination of the 
roughness exponent a. The following subsections are 
structured as follows. In A) we introduce the geo- 
metric elements or blocks (equivalent to the Kadanoff 
blocks in standard RSRG methods) suitable to deal with 
anisotropic situations. In B) we discuss the effective dy- 
namics of the previously defined blocks at a generic scale. 
In C) we introduce the RG equation and explain how the 
roughness exponent is determined. Finally in D) we an- 
alyze critically the approximations involved in general in 
the method. 



A. Geometric description 

The first nontrivial problem in the development of a 
RSRG approach is to find a sensible description of the 
geometry of the growing surface at a generic scale, i. e. 
how to build the analog of a block-spin transformation 
p3[ . Given the anisotropy of the system, the shape of 
the blocks must depend on the scale. Therefore, sub- 
dividing a cell in subcells is not a feasible task and the 
explicit construction of the block-spin transformation is 
not possible. 

Hence we develop an alternative strategy. To obtain a 
description at a generic scale k of the growing surface, we 
consider a partitioning of the (d-t-l)-dimensional space in 
cells of lateral size Lk = Lob'' and vertical size hk- Here 
6 is a constant and k labels the scale (Fig. |l]) . 

A cell is declared to be empty or filled according to a 
majority rule. In this way we pass from the microscopic 
description h{x, t) to a coarse-grained one at scale /c, fully 
defined by the number h{i,k,t) of filled blocks in the 
column i. Heights at scale k are measured in units of hk- 
The only characteristic vertical length at scale k is that 
fixed by typical intrinsic fluctuations of the surface of a 
lateral size Lk- This suggests to take 

hk = V~cW{Lk) ^ Lt, (5) 

where ^/c is a proportionality constant that will be dis- 
cussed later. This equation expresses the requirement 
of scale invariance in the geometric description. Any 
other choice would result either in a redundant descrip- 
tion (if hk/W{Lk) ^ as k oo) where too many (in- 
finite) blocks would be needed to describe fluctuations 
in the same column, or in a too coarse description (if 
hk/W{Lk) cx) as fc — > oo). By imposing Eq. (||), we 
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always have a meaningful covering of the surface upon 
scale changes. Observe that since in general a 7^ 1, the 
shape (i.e. the ratio of vertical to horizontal length) of the 
blocks changes with the scale k. Contrarily to the usual 
RG approach, the definition of the block-spin transfor- 
mation depends explicitly on the roughness exponent a, 
the calculation of which is the final goal of the method. 

The constant c in Eq. (^ fixes the unit of measure of 
our blocks. Its optimal value can be determined as fol- 
lows. The distribution of microscopic height fluctuations 
within a block can be mapped into an effective distri- 
bution with the same average h and standard deviation. 
For simplicity we take it to be bimodal 

P[h{x)] = p5{h{x) -[h+{l~ p)hk]} 

+ a-p)5{h{x)~[h-phk]}. (6) 

This distribution results from mapping all points with 
microscopic height larger than h to h+(l—p)hk and those 
smaller than h to h — phk- The parameter p describes the 
degree of asymmetry of the distribution: The fluctuations 
inside a block can then be calculated, using Eq. (|^), as 

W^{Lk)^p{l-p)hl (7) 

which implies that the constant c is given by 



For a symmetric distribution p —1/2 and therefore c = 4. 
In general the height distribution is not symmetric, i. 
e. there is some nonvanishing skewness and one must 
consider c 7^ 4. 

B. Dynamic description 

The second step in the construction of the RG pro- 
cedure is the definition of the effective dynamics at a 
generic scale k, i.e. the determination of the growth rules 
for the blocks defined in the previous subsection. The ef- 
fective dynamics will depend on a set of scale dependent 
parameters. The changing of scale induces a flow in the 
parameter space whose fixed points correspond to the 
scale invariant dynamics. 

Analogously to what happens in the usual application 
of the RG approach to equilibrium systems, it may hap- 
pen that mechanisms not appearing in the microscopic 
rule are generated upon coarse-graining. In the language 
of equilibrium systems this means that operators not in- 
cluded in the bare Hamiltonian can be generated itera- 
tively. Conversely, microscopic ingredients can prove to 
be irrelevant and be progressively eliminated when going 
to coarser scales. Therefore, exactly as in the equilibrium 
case the choice of the parametrization of the effective dy- 
namics is not trivial: Principles as the preservation of 
symmetries and conservation laws must be the guidelines. 



In general, the effective dynamics will be defined in terms 
of the transition rates for the addition of occupied blocks 
at a generic coarse-grained scale, that is, 

r[h{i,k)^h'{i,k)]=r{xlxl...,xl). (9) 

The number of parameters x\. is in principle arbitrary, 
although in the applications presented below it will be 
limited to one It is clear that the more complete 

the parametrization the better the final description of 
the statistical scale invariant state. We will discuss this 
problem in detail in subsection D. 

C. The RG Equations 

So far we have defined the geometrical and dynami- 
cal aspects of the coarse-graining procedure. These give 
us the necessary ingredients to introduce the RG trans- 
formation. The explicit derivation of it is based on the 
following property of the roughness W . Let us consider 
a d-dimensional system of linear size L and partition it 
in [L/hy^ blocks of size h'^ (labeled by the index j). It is 
straightforward to verify that the total roughness can be 
decomposed as 

+ (iT^P E m-h]\ (10) 

where h{j) is the average height within block j. The 
interpretation of this formula is simple: The first term on 
the right hand side is the averaged value of the roughness 
within blocks of size 6'', while the second term is the 
fluctuation of the average value of h among blocks. 

In our coarse-graining procedure this property is read 
as follows: If one takes L = Lk+i = hLk the first term on 
the right hand side is W'^{Lk), the total roughness within 
a block of size L^; the second is the roughness of the 
configuration in which blocks of size are considered 
as flat objects. This second contribution is obviously 
proportional to the square of the height of a block h\. 
Hence, employing Eq. (H), 

W\Lk+i) = W\Lk)+uj\b, k)hl 
= [l + CLo\b,k)\ W\Lk) 
= Fbik)W^{Lk), (11) 

where uj'^{b, k) is the roughness in the stationary state of 
a system of b"^ sites of unit height that evolves according 
to the dynamical rules specified by (x^, x|, . . . , a;^), and 

Fhik)=[l + cu;^ib,k)]. (12) 

Note that the dependence on the scale k is only through 
the parameters {xk}- 
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Eq. (|7Jj is the equation that relates the width at scales 
k and k+1. In order to proceed further, we must evaluate 
the function Fb{k), or equivalently a;^(6, k). To do so, we 
identify all the possible surface configurations of a system 
of composed of sites, and write down a master equation 
for their associated probabilities pi 



(13) 



Pi^j is the rate for the transition between configuration 
i and j and depends on the set of parameters {xk\- Im- 
posing the stationarity condition dtPi = and the nor- 
malization Pi = 1 the master equation can be solved. 
If we call Wf the roughness of configuration i, then we 
can write 



(14) 



Depending on the particular structure of the master 
equation the explicit solution of the previous equation 
may be difficult or impossible. In such cases it may be 
more useful to determine uj'^{b,k) numerically by per- 
forming (relatively small) Monte Carlo simulations. We 
will describe examples of both analytical and numerical 
computations of uj'^{b, k). 

Let us suppose now that uj'^ (6, k) has been determined. 
Eq. (^l|) gives an explicit relation between the roughness 
at two different scales. Observe that so far the scale in- 
variance idea has not been implemented. We have just 
studied how the width changes upon changing the level of 
description. The last task to be performed is the determi- 
nation of the RG transformation relating the parameters 
of the dynamics at scale k with those at scale k + 1. This 
is done by means of a self-consistency requirement for 
the description of the same system at two different lev- 
els of detail, i.e. the total width of a system should be 
independent on the size of the blocks we use to describe 
it. To make this idea more precise, let us consider the 
case of a dynamics parametrized by only one parameter 
Xk- Let us take a system of size L = Lk+2- By applying 
Eq. (|ll|) we have 



W^iLk+2) = Fbixk+i)W^iLk+i). 



(15) 



This procedure can be iterated again on each of the re- 
sulting systems of size Lk+i, obtaining 



(16) 



The same quantity can alternatively be computed by con- 
sidering directly the whole system as composed by b^^ 
systems of size Lk- Applying again Eq. (O) 



W^iLk+2)^F,2{xk)W^iLk). 



(17) 



Imposing the consistency of the two procedures one has 
an implicit RG transformation for Xk 



Fb{xk+i) = 



Fb^jxk) 
Fb{xk) ' 



or explicitly 



Xk+i = R(xk) = F^^ 



Fb2{xk) 



Fb{xk) 



(18) 



(19) 



This equation provides the evolution of the parameter 
under a change of scale. 

If a fixed point x* such that 



R{x*) 



(20) 



exists, then the parameter x* characterizes the scale in- 
variant dynamics of the system. The knowledge of it di- 
rectly allows the determination of the exponent a. Since 
W'^{Lk+i)/W^{Lk) is equal to 6^" we have 



a = lim i logf, 



lim \ \o%bFb(k) = \ \o%bFb{x*). 

k—*oo A A 



(21) 



To analyze the stability of the fixed point we linearize 
the RG transformation around it 



Xk^\ 



Rixk) 



R'{x*){xk-x*). (22) 



Hence if |i?'(x*)| < 1 the scale invariant dynamics spec- 
ified by X* is an attractive fixed point under changes of 
scale. 

Extension of the previous formalism to the case of n 
parameters of the dynamics is straightforward. The n 
RG transformations are obtained by imposing the consis- 
tency of the description of the same system when divided 
in 2'* and 4'' blocks, in 4"* and 16*^ blocks, and so on [n7|. 



D. Approximations 

Let us discuss now the approximations involved in the 
method. There are two steps where approximations come 
into play: The first is the choice of the parametrization 
of the scale invariant dynamics. The second is the com- 
putation of uJ^{b, k). 

With respect to the first problem, it is reasonable to 
expect that under coarse-graining the microscopic dy- 
namics will flow towards a scale invariant dynamics de- 
pending in principle on an infinite number of parameters. 
This proliferation is analogous to what happens in RSRG 
approaches to equilibrium systems. The restriction to a 
finite (and small) number of parameters involves unavoid- 
ably an approximation, due to the projection of the RG 
flow onto the sub-space spanned by these parameters. 

However, a very important difference with respect 
to equilibrium critical phenomena is that here the 
scale invariant dynamics is "self-organized critical" , that 
is, there are no relevant operators. Only irrelevant 
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fields, with negative scaling dimensions, need to be 
parametrized. The system is by definition on the crit- 
ical manifold and, by iteration of the RSRG transforma- 
tion, it converges to the stable fixed point, without any 
fine tuning of parameters. The projection onto a low- 
dimensional parameter space yields a projected RG flow 
which will share these same properties. The fixed point 
in this sub-space, being the projection of the actual fixed 
point in the high-dimensional space, will have the same 
qualitative properties. Even the simplest parametriza- 
tion capturing the correct symmetries of the dynamics 
can provide a quite accurate determination of the prop- 
erties of the system in this case. On the contrary, when 
relevant fields are present, as in second order phase tran- 
sitions, truncation effects are quite dramatic. The reason 
is that relevant fields have, in general, a non-vanishing 
component on any discrete (lattice) operator. Any ap- 
proximation due to truncation is amplified by the RG 
iteration thus driving the flow out of the critical mani- 
fold along the relevant directions. The determination of 
the fixed point becomes then very difficult [Q. 

The second source of approximation is the computa- 
tion of k). As stated above this quantity is the sta- 
tionary roughness of a system composed of b'^ substrate 
sites evolving according to the dynamical rules specified 



by the parameters {x\, . 



This is a perfectly well 



defined quantity that may in principle be computed to 
any degree of accuracy by solving the master equation. 
However very often the structure of the master equation 
is too complicated to allow for a full solution. One then 
has to devise suitable simplifications to make the analyt- 
ical computation feasible. This involves approximations 
that affect the final result. We will see an example of 
this way of proceeding and discuss how the effect of the 
approximation can be controlled. Alternatively, when b 
and d are not too large, one can resort to the numer- 
ical evaluation of uj'^{b,k). In practice this boils down 
to performing Monte Carlo simulations of small systems 
evolving with different values of the parameters {2;^}. It 
is important to stress that the MC procedure involves 
no approximation, except for the fluctuations associated 
with statistical sampling. We will describe below an ex- 
ample of this alternative way of computing w^(6, k). 

A delicate issue is also the choice of the boundary con- 
ditions. In the conceptual framework described above, 
uP{b^k) is the roughness of a section of size b of an in- 
finitely extended surface. This would suggest the use of 
open boundary conditions. On the other hand, when in- 
tegrating out degrees of freedom relative to height fluctu- 
ations inside the cell, one should not consider the fluctua- 
tion of the average slope. This slope effect is eliminated if 
one uses closed (i.e. periodic) boundary conditions. Even 
though the choice of the appropriate boundary conditions 
is not trivial, we will see, in the KPZ case, that the use of 
periodic or open boundary conditions has little effect on 
the value of the exponent. Furthermore, one expects that 
both truncation errors and those induced by neglecting 
fluctuations of boundary conditions vanish as the param- 



eter b grows. Arguments in support of this conclusion are 
reported in the Appendix I. 



III. RG FOR KPZ DYNAMICS 

A. The problem of KPZ growth 

The Kardar-Parisi- Zhang equation is the minimal con- 
tinuum equation capturing the physics of rough surfaces. 
After appearing in 1986 an overwhelming number of 
studies has been devoted to elucidate its properties [^,|| . 
It reads [il 



dh{x, t) 
dt 



■q{x,t). 



(23) 



where /i(x, t) is a height variable at time t and position x 
in a c?-dimensional substrate of linear size L. v and A are 
constants and 77 is a gaussian white noise. As a conse- 
quence of a tilting (Gafilean) invariance a-l-z = 2, 
and since in general z = a/ (3, there is only one indepen- 
dent exponent, say a. The difference between the KPZ 
equation and the linear equation (Edwards- Wilkinson), 
describing surfaces growing under the effect of random 
deposition and surface tension, is the presence of a non- 
linear term proportional to A. This nonlinear term is 
generated by microscopical processes giving rise to lat- 
eral growth, i.e. the fact that growth velocity is normal 
to the local surface orientation. 

Exact results |^J^ indicate that in d = 1 there is only 
a rough phase for KPZ with a — 1/2. Instead, standard 
field theoretical methods predict the presence of a rough- 
ening transition above d = 2 [ p^ ; i.e., there are two RG 
attractive fixed points and an unstable fixed point sepa- 
rating them. More specifically, there is a gaussian fixed 
point with a = describing a flat phase (characterized 
by a vanishing renormalized nonlinear coupling) and a 
nontrivial one describing the rough phase (in which the 
renormalized nonlinear coupling diverges in perturbation 
theory). Perturbative methods fail to give any predic- 
tion for the exponents in the rough phase. For d > 2, 
an e-expansion (d = 2 -I- e) around the gaussian solution 
can be performed and the exponents at the roughening 
transition evaluated to all orders in perturbation theory 
|p0|,p| . These results seem to indicate the presence of an 
anomaly in d = 4 for the roughening transition. This has 
been interpreted as an indication that dc = 4 is the upper 
critical dimension for the rough phase, i.e., for the strong 
coupling fixed point [pT[-p3|. Above this dimension the 
exponents should take the values known for d = 00 [ p^ . 
Applications of non-perturbative methods such as func- 
tional renormalization group and Flory-type argu- 
ments also suggested that dc = 4, in agreement with 
a 1/d-expansion [ p6| around the d = 00 limit. The mode- 
coupling approximation led to contradictory results, sug - 
gesting the existence of a finite dc or dc = 00 [p7[. 
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Arguments for a finite dc based on directed or inva- 
sion percolation have also been proposed. 

On the other hand, numerical results seem to indicate 
that the exponent a decays continuously with the system 
dimensionality up to d = 7, excluding therefore d — A as 
upper critical dimension [ pO| . 

Finally, some doubts have been cast on the validity 
of the continuum approach to study rough surfaces [pl| . 
Summing up, the issue of the behavior of the KPZ dy- 
namics for d > 2 is a highly debated one, and it is ex- 
tremely desirable to have alternative approaches shed- 
ding light into the problem. In what follows we present 
the application of our new RG scheme to KPZ growth. 



B. Simplest RG scheme 

1. Parametrization of the dynamics 

The modelization of the dynamics at a generic scale 
should keep the number of parameters to a minimum and 
catch all the relevant physical mechanisms of the process. 
The main feature of the KPZ dynamics is lateral growth. 
Therefore we take as the only parameter defining the dy- 
namics at a generic scale fc, the ratio Xk between lateral 
and vertical growth (i.e. random deposition). More for- 
mally, the growth rate for the addition of an occupied 
block on column i is 



r[h{i)^h{i) + l\ 

1 -I- Xk ^ max [0, h(j) — h{i) 



(24) 
(25) 



jn.n.i 



Eq. ( p5| ) states that the rate for lateral growth is pro- 
portional to the difference in height between neighboring 
columns (Figj2). Overhangs, known to be irrelevant on 
large scales ]16[ |, are not allowed. This dynamics can be 
seen as a generalization of the Eden growth model. 

Few observations are in order. We call the parameter 
X appearing in Eq. (25) "lateral growth" parameter, but 
this is an abuse of language: x^ cannot be identified with 
the parameter A of the KPZ equation. Instead the term 
that multiplies Xk in Eq. ( p5| ) is a combination of the dis- 
cretized Laplacian, of the discretized square gradient and 
of other discrete operators. The explicit dependence of x 
on V and A cannot be disentangled. Other parametriza- 
tions are clearly possible and they will be discussed be- 
low. 

Eq. (^) has the nice feature that it contains as limiting 
cases both the random deposition process (x^ = 0) and 
the infinitely strong "lateral growth" [x^ = oo) leading to 
fiat surfaces. Most importantly, it is easy to see that x* = 
oo is, by construction, a fixed point of the RSRG with 
a = 0. This feature makes it possible the determination 
of the upper critical dimension above which the stable 
solution leads to a = 0. In this situation we expect 
x* = cx) to be an attractive fixed point. Below the critical 



dimension, on the other hand, the fixed point x* = oo 
must be unstable and an intermediate fixed point with 
finite a must appear. The RSRG accommodating for a 
fixed point at a;* = oo, naturally allows to address the 
issue of the upper critical dimensionality. 



We restrict ourselves for the moment to the one- 
dimensional case and illustrate in detail the application 
of the RG approach, i. e. the computation of (6, fc), the 
determination of the scale invariant dynamics and of the 
exponent a. It is very instructive to consider first the dy- 
namics Eq. (I25I ) supplemented by the condition that the 
height difference between adjacent columns is restricted 
to values such that (|A/i| < Ahmax), with Ahmax = 1- 
This greatly reduces the number of possible surface con- 
figurations allowing for a full analytical treatment. For 
the system of size b — 2, assuming periodic boundary 
conditions, there are only 2 nonequivalcnt configurations, 
while for the system of size 4 there are 6 of them (Fig. ||). 
Using the definition Eq. ( p5|) of the growth rates for the 
addition of a block, one has simply, for 6 = 2, 

Pi^i = 
Pi^2 = 2 
P2^1 = l + 2xk 

P2^2 = 0. (26) 

In configuration 1 only vertical growth is possible (in two 
sites) leading always to configuration 2. Only one site 
can instead grow in configuration 2 and the rate for this 
is the sum of the rate of one vertical and two lateral 
contributions. Hence the master equation reads 



^ = (1 + 2xk)p2 - 2pi 
ot 

^ = 2pi - (1 + 2xk)p2- 
ot 



(27) 
(28) 



Imposing the stationarity condition dtpi — dtP2 — and 
the normalization pi + /)2 = 1 one has 



P2 



2, + 2xk 



(29) 



and then considering that the width associated with con- 
figurations 1 and 2 is, respectively, and 1/4 



L0''{2,Xk) 



4(3 + 2xk) 



(30) 



For the system of size 6^ = 4 one finds in an analogous 
way 



w2(4,Xfe) 



51 + mxk + 404 

4(47 + 106a;fe + 68a;2+ 8x1) ■ 



(31) 
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Plugging these two expressions into the RG equation jl^ ) 
with = 40 one finds that the exphcit form of the RG 
transformation is 



R{x) = 



293 + 804.x + 636x^ + 160x'^ + 32a;^ 



2(59 + 148j: + 156a;2 + 64x3) 
and that there exists only one finite fixed point for 

X* ~ 2.08779 . . . 
Such a fixed point is attractive since 
R'{x*) ~ -0.03548... 



(32) 

(33) 
(34) 



Hence no matter how small or large the microscopic value 
of X is, upon coarse-graining the dynamics flows towards 
an attractive scale invariant dynamics, characterized by 
a ratio x* of the lateral to vertical growth rates. The 
roughness associated with this scale invariant dynamics 



summarized graphically in Fig. |[ We find a finite at- 
tractive fixed point for all dimensions, with an exponent 
a in remarkably good agreement with the best numerical 
results available This is the first theoretical ap- 

proach providing estimates for the roughness exponent 
that match in all dimensions with numerics. No anoma- 
lies are found for d = 4 where other approaches find an 
upper critical dimension. The extrapolation to d — > oo 
suggests that the fixed point is always stable and that a 
decreases with d but remains always nonvanishing. The 
fixed point parameter x* grows exponentially with the 
dimension. These results are confirmed by an analyti- 
cal expansion of the method in high dimensions, that is 
presented in subsection D. 



C. Robustness of the results 



1. b>2 



a = -log2i^b^2(a;*) ^ 0.177352. 



(35) 



that must be compared with the known exact value a = 
1/2. 

The apparent poor performance of the method is due 
to the assumption that Ahmax = 1, which allows for full 
analytical treatment, but is clearly wrong. The point 
is that even if at the microscopic level the dynamics 
is of restricted type, the effective dynamics at generic 
scale defined by the renormalization procedure will pro- 
liferate in a non restricted one. Allowing larger steps 
{Ahmax > 1) increases the number of superficial config- 
urations and makes the analytical determination of the 
function uj'^{b,Xk) impossible. Still this task can be per- 
formed numerically via simulation of systems of such a 
small size. Fig. || reports the results obtained by consider- 
ing increasing values of Ahmax- The value of x* (Ahmax) 
converges already for Ah„iax = 8 to x* ~ 0.726, corre- 
sponding to a value a = 0.507 ... in excellent agreement 
with the exact value. Further increases of Ahmax do not 
change the results, indicating that in the scale invariant 
dynamics the probability of steps larger than 8 is negli- 
gible. 



3. d>l 

The computation via Monte Carlo method of cli^(2, Xk) 
and w^(4,Xfc) can be performed with very little compu- 
tational effort also in higher dimension. We considered 
d = 1, ... ,9 with less than a week of CPU time of a 
workstation. The results are reported in Table I and 



In order to analyze the stability of the results upon 
increasing the value of b it is convenient to introduce the 
function 



(36) 



With this definition one can express the fixed point con- 
dition Eq. (M) as 



ab{x*) = Q!f,2(a;*), 
and see that the fixed point is stable if 



\R{x* 



a'f,{x*) 



< 1, 



< 



V^{x*) 



< 1. 



(37) 



(38) 



(39) 



^In d = 1 the distribution is known to be symmetric. 



Such a formula can be extended also to the case where 
the size of the larger system considered is not 6^ but a 
generic b' > b. We study the stability of the results for 
growing b by computing ab{x) with bi = 2,4,8,16,... 
and imposing the consistency between two successive bi. 
The value of b indicated in the plots is the smaller one; 
for instance, 6 = 4 label the results obtained imposing 
the consistency between 6 = 4 and 6' = 8. 

In Fig. ^ we report the plot of the curves ab{x) in d = 1 
for 6 = 2,4,8,16,32. Remarkably they all meet practi- 
cally at the same point indicating that x* and a virtually 
do not change by increasing the number of cells. Fig. |^ 
reports the values of the exponent a in d = 1 (empty 
circles). Observe that fluctuations are extremely small. 
The value of the fixed point parameter x* is reported 
in Fig. P (empty circles). Again it remains practically 
unchanged when 6 grows. 
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In higher dimensions the results are less stable. The 
values of a and of x* for d = 2, 3 and 4 are reported in 
Fig. and ^ respectively (empty symbols) . A clear trend 
is present for d — 2: the exponent initially decreases as b 
is increased, then reaches a minimum and starts growing. 
This behavior is reflected in the value of x* that first 
grows and then decreases. 

The decreasing part of the pattern is present in the 
analogous plots for d = 3 and d — A. For large dimen- 
sions however, it is increasingly more time consuming to 
perform the computation for large systems. In particular 
for d — 4, the largest system that could be simulated is 
6 = 16 and for such a system size the trend is still de- 
creasing. Therefore it is not possible to decide from a 
numerical point of view whether for larger values of b a 
would converge to zero or to a finite value. These data do 
not provide any conclusive indication on whether d = 4 is 
the upper critical dimension for KPZ growth. However, 
as it will be shown below, such a conclusion is ruled out 
by the results with other parametrizations and by the 
analytical large-d expansion of the method. 

The reason for the difference in the stability of the 
results for large number of cells in d = 1 and d > 2 
is probably related to crossover phenomena. In the RG 
flow there are two competing fixed points; this reflects the 
existence of two universality classes, strong coupling KPZ 
and EW. In d = 1 this fixed points are associated with the 
same roughness exponent a = 1/2 and to similar values 
of X* . Therefore, any crossover phenomenon between the 
two fixed points has little effect in our formalism. In d > 
2 instead, the two scale invariant dynamics are associated 
with different exponents and also very different values of 
the parameter x* , which is finite for KPZ and infinite for 
EW. We interpret the initial decrease in the value of a 
in the KPZ case as the effect of a crossover caused by 
the presence of the EW fixed point. It is not clear to us, 
however, why the fixed point found for 5 = 2 is so close 
to the results of the numerical simulations. 



2. Open boundary conditions 

The calculation of uj'^{b, k) can also be performed with 
open boundary conditions, that is assuming that the 
height of the columns outside the system which are in 
contact with the boundary is the same of their neighbors 
inside the system. This means that no "lateral growth" 
event can be caused in the system by the environment 
around it. 

The results are also presented in Fig. ^ and ^ (filled 
symbols). Interestingly, in this case the accuracy of the 
method for 6 = 2 is not as good as for periodic boundary 
conditions, but the error remains below 10%, indicating a 
low sensitivity to the boundary conditions even for small 
number of cells. For larger number of cells the difference 
goes quickly to zero. 

For higher dimensions the general dependence of a;* 



and a on 6 remains unchanged: In d = 2 a is initially 
high, then decreases and finally increases again. The vari- 
ations with b are however less strong than when periodic 
boundary conditions are considered. For d = 4 it is more 
clear than in the case with periodic boundary conditions 
that a does not converge to zero for large b. 



3. Other parametrizations of the dynamics 

As stated above the parametrization ( p5| ) of the KPZ 
scale invariant dynamics is by no means unique. Actually, 
given the problems of slow and nonmonotonic conver- 
gence towards the asymptotic values, it turns out clearly 
that the parametrization (|2^) is quite far from being opti- 
mal and better parametrizations would help. In order to 
keep things as simple as possible we started considering 
transition rates of the form 



(40) 



jn.n.i 



with 7 constant; for 7 = 1 it coincides with (25). By 
comparing the values of a obtained with several 7 an in- 
teresting pattern can be spotted (Fig. ^ . While for small 
number of cells b the estimate gets worse with increasing 
7, the opposite is true for large b. For large values of 
7 the estimate for a converges quite rapidly. For 7 = 9 
and 7 = 20 we find on the largest systems a — 0.399, 
suggestive of a convergence towards 0.4. For d = 4 the 
sizes that can be simulated are too small to allow the de- 
termination of the asymptotic value of a. However, it is 
clearly seen that a does not go to zero as b is increased. 

The same type of behavior is found by using an expo- 
nential parametrization of the dynamics 



r,: = 1 



Xk 



J2 exp{7.max[0,Mj)-M*)]}-l- (41) 



In d = 2 for large 7 the estimate of a on the largest 
system is 0.399, exactly as with Eq. (^o|). In d = 4 again 
we cannot precisely determine where a is converging to. 
Again the data strongly suggest that this limit is finite. 

The study of these two alternative parametrizations of 
the dynamics consistently indicates a value of a = 0.399 
in d = 2 and a finite a > in d = 4 suggesting strongly 
that 4 is not the upper critical dimension of the KPZ. 

One could in principle imagine a parametrization of 
the effective dynamics, more in the spirit of the KPZ 
equation, of the type 



ri = 1 + iyk\V^h{i^ \ + \k[Vh{i)f 



(42) 



However, there is no reason for believing that such a 
parametrization would be better for the KPZ rough 
phase; additional operators are very likely to be present 
in the scale invariant dynamics. Moreover the dynamics 
described by Eq. ( ^ ) is plagued by numerical instabili- 
ties, as pointed out by Bray and Newman |32|. 



8 



D. The d ^ oo limit and the upper critical dimension. 

The results presented so far show that a > even for 
large number of cells in d = 4, thus indicating that 4 
is not the upper critical dimension for KPZ growth pro- 
cesses. By using the RG procedure it is actually possible 
to go beyond this numerical conclusion: The existence of 
any finite upper critical dimension can be ruled out. This 
result is obtained when the function uj'^{b,Xk) is com- 
puted analytically in the large-d limit. The basic fact 
allowing this calculation is that when d^ 1 one expects 
a <C 1, which suggests that surface fluctuations are small 



uj{b,Xk) - 6" ~ l + alnb + Oia^). 



(43) 



For small b one may reasonably account for the fluctu- 
ations of the interface by considering only two possible 
values of h{i), ho ("low sites") and ho + 1 ("high sites") 
(Fig. |ll|). Starting from a flat surface {h{i) — 0,Vi), one 
considers growth events occurring according to the rates 
Eq. (p5|), with the restriction that no block can be de- 
posited on top of an already grown one. Only when the 
whole layer at height 1 is grown one allows growth to 
level 2 and so on. This approximation allows the analyt- 
ical evaluation of oj^{b, k), the identification of the fixed 
points and the study of their stability. We will check a 
posteriori the consistence of the results with the assump- 
tion, and see that the existence of a finite upper critical 
dimension can be excluded. Let us now present the de- 
tails of the calculation. 

Within the "two layers" approximation it is convenient 
to group together all configurations with the same num- 
ber of high sites: we will call "state" n the set of all 
surface configurations with n sites at height ho + I and 
the remaining b'^ — n sites at height ho- The state n = 0, 
corresponding to a flat surface, is equivalent to the state 
with n — b'^. This classification is useful because the only 
transitions permitted from state n are those to state n~\-\. 
The master equation for the probability pn of being in 
state n (i.e. of having any of the configurations with n 
high sites) is then greatly simplified 

dtPn = Pn-ir{n - 1 n) - p„r(n -> n + 1). (44) 

r{n ^ n -I- 1) is the average of all the rates ( p5| ) for the 
growth processes that transform one configuration with 
n high sites in one with n + 1 of them. We can write this 
quantity in the form 



r{n 



+ 1) = [b'^ - n) + Xk^r,. 



(45) 



The first term on the right hand side is simply the total 
rate of vertical growth [1 in Eq. (|25|)] for configurations 
with n high sites. Observe that it is obviously equal to 
the number (6'^ — n) of sites where vertical growth is al- 
lowed. Xk^n is the rate for lateral growth: f2„ is the 
average number of lateral walls in configurations with n 
high sites. Its precise computation is not easy, since it 



would require the knowledge of the stationary probabil- 
ity for each configuration belonging to state n. However, 
when Xk — the computation is trivial since growth oc- 
curs only via uncorrelated deposition and high and low 
sites are randomly distributed. The number of low sites, 
where growth is allowed, is b"^ — n; each of them has 
2d neighbors which are occupied with probability n/b'^. 
Hence the average number of lateral walls is 



(b'^ - n)2d 



(46) 



The form of Qn for a; 7^ is in general more compli- 
cated, but a numerical computation for large dimensions, 
namely for d = 7, shows that Eq. ( fig) is a good approx- 
imation for all values of Xk (Fig. ]l2|). We assume the 
validity of Eq. (E^) for all values of x^- This leads to 



r{n ^ n + 1) = {b'^ - n) 

The stationary solution of Eq. (| 
r(0 -> 1) 



l + 2d^x. 



Po 



r{n 



1)' 



(47) 



(48) 



n=0 



Pn 



By imposing the normalization condition ^ 
and approximating the sum by an integral one obtains 

Ld r / , , „ , X -1 ^ -1 



1 



Po= 1 



2dxk 



2d In & + In 



1 + 2dxk 
b'^ + 2dxk 



(49) 



Equations (^) and (^9|) provide a complete descrip- 
tion of the stationary probability density. Given that 
the roughness of all configurations with n high sites is 
n/b'^{l — n/b'^), the total roughness of the surface can be 
computed as 



u: {b,Xk)^ }_^Pn 



(50) 



Using the fact that b'^ ^ 1 and assuming dxk » 1, we 
obtain 



uj^{b,Xk) = Po 



2dxk 



(51) 



Inserting Eq. ( pl| ) with 6 = 2 and 6 = 4 in the fixed point 
equation (pOl) yields, to leading order in d, 



■)d+l 



In 2. 



(52) 



The assumption dxk 3> 1 is therefore self-consistent for 
sufficiently large d. Notice that an exponential depen- 
dence of a;* on d was already found in the numerical im- 
plementation of the method (Fig. ||). Using Eq. ( pl| ) one 
obtains the value of the roughness exponent 



1 1 



3(ln2)2 d' 



(53) 
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Finally, by computing the derivative of the RG transfor- 
mation at the fixed point 



R'ix*) = -1 



1 



21n2d 



0(l/d'), 



(54) 



we see that the fixed point is attractive for all finite di- 
mensions. 

In conclusion we find that for large-c? the RG has a 
fixed point x* corresponding to an exponent a ~ 1/d and 
therefore strictly positive in all finite dimensions. On the 
contrary the existence of a finite upper critical dimension 
would have implied, for d > dc, either the absence of a 
finite fixed point or its instability. 

At this point we must use the analytical result to check 
the consistency of the two layers assumption. Such as- 
sumption is correct provided the rate of its violation is 
negligible for all values of n. Processes violating the as- 
sumption are those in which an event of vertical growth 
occurs on top of an high site. Their rate in state n is 
^up — n, that must be compared to the total rate of 
processes not violating the restriction r(n — > n + 1) com- 
puted for X — X* . By imposing r„p(n) ^ r{n ~^ n + 1) 
we get 



n < (fo'' - n) 1 + 



2dn 



(55) 



Let us consider n = — 1 which is the situation that 



maximizes r^p and minimizes r{n 



1). Then 



6'' - K 1 



2d{¥ - 1) 



Since 6"^ ^ 1, this means 

b'^ < 2dx* - 2'^+'^d. 



(56) 



(57) 



Hence the two layers assumption is correct for b — 2 
but fails for 6 = 4. Therefore the value of uP'iA^Xk) is 
systematically underestimated by Eq. ( ^l| ) since fluctu- 
ations involving more than two layers are neglected de- 
spite being likely. The consequence of this on our results 
is understood by considering Eq. (|8|) . In such a formula 
we estimate correctly the left hand side, while the right 
hand side is underestimated (Fig. Since the fixed 

point parameter x* and the exponent a are given by the 
intersection of the curves it is clear that we get an upper 
bound for x* and a lower bound for a. This is confirmed 
by the comparison of our estimate of a, Eq. (|5^), and 
X* with the numerical results of Ala-Nissila et al. and 
the value of x* computed numerically for d — 1, . . . , 9 
(Fig. Hi). 

These results have been obtained for the smallest val- 
ues of 6, namely 5 = 2. In the previous sections we 
showed that for low dimensions the results for small b are 
in good agreement with numerics, but for larger b there 
are deviations. A very reasonable question therefore con- 
cerns the robustness of the large-d results when b grows. 



As we have shown above the two layers approximation 
breaks down for b > 2. In order to extend the above 
calculation to larger b one should replace the two layers 
approximation with some less restrictive but still doable 
calculational scheme. We have not been able to fulfill this 
task and hence we cannot directly show whether a fixed 
point exists for finite x* when b — > oo. However, the two 
layers assumption is valid for any b in the neighborhood 
of the fixed point x* — oo that gives a flat surface a = 0. 
This fixed point exists in any dimension, and its stabil- 
ity can be safely analyzed using the previous two layer 
assumption as follows. 

Let us introduce e = 1/x. The derivative of the RG 
transformation at the fixed point e* = is (see Eq. (|3^)) 



R'(e = 0) 



K(e - 0) 

(6 = 0) 



1 



(58) 



where now prime indicates derivative with respect to e. 
To first order in e we have (see Appendix II) 



Then 



oU. = 0) = ^c,.--. 



Analogously 



a',{e = 0) 



2 In 6 



d+l 



Hence 



R'{e = 0) = 6' 



d+i 



1 > 1 



(59) 



(60) 



(61) 



(62) 



and the fixed point corresponding to a = is unstable. 
As a consequence, a finite fixed point with a > must 
exist and be stable when 6 — > cx) for any large and finite 
d. This supports the conclusion that there is no finite 
upper critical dimension. 



IV. RG FOR THE EDWARDS- WILKINSON 
DYNAMICS 

So far we have applied the new RG method to a KPZ- 
like dynamics. Now we intend to show that it is more gen- 
eral and can be applied to growth mechanisms belonging 
to universality classes other than KPZ. In particular we 
study in this section its application to the exactly solv- 
able, Edwards- Wilkinson (EW) equation for which the 
roughness exponent is known in any dimension. In par- 
ticular a — 1/2 in d = 1, while a = for > 2 with 
logarithmic corrections at d = 2. 

The parametrization Eq. ( ]2^ ) of the dynamics de- 
scribes a growth model where only deposition events can 
take place and the symmetry between up and down in 
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the h direction is clearly broken. Such a dynamics is 
inherently out of equilibrium and therefore cannot ac- 
commodate the scale invariant dynamics of the Edwards- 
Wilkinson growth process, which is an equilibrium one, 
with growth rules symmetric along the growth direction. 
We now introduce a generalized dynamics which admits 
the KPZ and EW dynamics as particular limiting cases. 
Let us consider the quantities 

Kail) = ^ max[0, - 

jnni 

Ku{i) = max[0, h{i) - h{j)]. (63) 

jnni 

In the KPZ case described so far we have allowed only 
deposition of particles and written 



and as J^i Kuii) = J2i ^d{i), we have that 



(64) 



We now allow also for evaporation of particles. That 
is, we consider the transition rate for site i as 

n^l + Xk\eKu{i)-il-e)Kd{i)\ (65) 
and with probability 

Pb = (66) 

a random deposition/evaporation event takes place {hi 
hi + 1 with probability e and hi—>hi — l with probability 
1 — e), while with probability 



(67) 



(68) 



Pi = Xk\eKuii) - (l - e)Kdit)\/n 
we have a "lateral" event 

, 6A-„(z) - (1 - e)Xrf(z) 
\eK4t)-{l-e)K4i)\- 

For e = 1 only deposition is allowed and we have the 
transition rates for the KPZ dynamics. For e = 1/2, we 
have up-down (deposition-evaporation) symmetry and 
the rates are 



n ^i + Xk\V^h{i)\ 



(69) 



where \/^h{i) = [Ku{i) — Kd{i)]/2 is the discretized 
Laplacian evaluated at site i. Therefore we expect this 
case to correspond to EW dynamics. Let us verify that 
for e = 1/2 the average interface velocity does not depend 
on the interface configuration (which is a basic property 
of EW dynamics) 



Id 



Since 



\Ku[i) - Ka(i)\ 

n \ \Ku(i) - Kd[i)\ 

- ix[Ku{i) - K^m 



Ah^ = Pb-0 + Pi- 



1 



(72) 



(71) 



With this generic dynamics we can perform the RG 
procedure exactly as in the KPZ case. The evaluation 
of the function uj'^{b,x) is carried out again using small 
Monte Carlo simulations with periodic boundary condi- 
tions. The results are reported in Fig. [ij. For d = 1 the 
value of a for & = 2 is ~ 0.4 below the exact value 1/2, 
but for b > 2 the correct value is rapidly approached. 
The situation is completely different in d = 2. In such 
case for 6 = 2 the exponent a is around 0.25, but when b 
is increased, the fixed point is shifted monotonically to- 
wards oo. In c? = 3 the behavior of the fixed point for 
finite x* is similar. From these plots we can conclude 
that the behavior of the EW dynamics is very different 
in c? = 1 and d > 1. For d — 1 there is a stable fixed point 
with a = 1/2. For d > I there is no stable fixed point 
with a 0. Hence the RG method is able to capture 
the difference between KPZ and EW dynamics and cor- 
rectly describes both the rough and the flat phases of the 
Edwards- Wilkinson growth. We speculate that the rea- 
son why one needs to consider 6 > 2 is probably related 
to the fact that on a system of size b = 2 the discretized 
Laplacian and square gradient take the same form. 



V. DISCUSSION 

In the previous sections we have introduced a general 
method for studying surface growth models by means 
of a real space renormalization group procedure. The 
anisotropy of the scale invariant properties of surfaces 
makes the very definition of the RSRG highly nontrivial, 
since the direct integration of degrees of freedom at small 
scale cannot be performed. For this reason we had to 
devise an alternative route: the main ingredient is that 
the integration of degrees of freedom is performed im- 
plicitly by imposing the self-consistency between two de- 
scriptions of the same object at different levels of coarse- 
graining. The application of such an approach to the 
KPZ dynamics yields several results that can be summa- 
rized as follows. 

• The scale invariant dynamics is identified and 
parametrized as a function of the "lateral growth" 
parameter x. This parameter turns out to have a 
non-trivial attractive fixed point under RG trans- 
formation for all dimensions. 

• The KPZ roughness exponent a estimated for small 
b is in very good agreement with large scale simu- 
lations of discrete models. 

• For larger values of b the estimate of a is stable 
in d = 1 while it changes noticeably for d > 2, 
presumably owing to crossover effects. When b 
oo it converges towards the correct result. 
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• The results are robust with respect to changes in 
the parametrization of the dynamics and in the 
boundary conditions. 

• No evidence is found of the existence of an upper 
critical dimension for KPZ. Moreover, we show very 
strong evidence that no such an upper critical di- 
mension exists. 

• By changing the nature of the parametrization of 
the dynamics at generic scale, the method is able 
to describe the EW dynamics and capture the ex- 
istence for it of an upper critical dimension above 
which only a trivial (flat) phase exists. 

Regarding the general nature of the approach it is 
worth remarking that the key point in the method is the 
identification of the scale invariant dynamics. In some 
sense the procedure can be seen as a kind of finite size 
scaling approach allowing for the evaluation of scaling 
exponents via the extrapolation of small size MC simula- 
tions. However the crucial point is that the MC data do 
not directly determine the exponent; they rather allow 
the identification of the scale invariant parameters of the 
dynamics which in turn determines the exponent. 

With respect to the estimates of the roughness expo- 
nent for small b, it is remarkable that the accuracy (in the 
sense of the discrepancy with known numerical results) 
seems to be the same in all dimensions d > 2. This is 
not the typical situation in ordinary critical phenomena, 
where usually RSRG methods fail in high dimensions. 
There are at least two reasons why usual RSRG schemes 
are inaccurate in high dimensions: 

a) The necessity of defining an explicit geometrical 
mapping between degrees of freedom at two scales (span- 
ning rule, majority rule, bond-moving, etc.). 

b) The presence of relevant fields and the need to com- 
pute exponents from the derivatives of the RG transfor- 
mation at the fixed point. 

Due to the success of field theory in high dimensions 
(close to dc) in usual critical phenomena, RSRG meth- 
ods have been mostly devised to work in low dimen- 
sions where the predictions of the e-expansion become 
less reliable. RSRG methods based on an explicit ge- 
ometric mapping (block-spin transformation) are quite 
accurate (and sometimes even exact) close to the lower 
critical dimension. Problems related with this geometric 
transformations become worse and worse as the dimen- 
sion increases. In our perspective, the only limitation 
has to do with the quality of the parametrization of the 
RG transformation. For example, the parametrization 
of the RSRG transformation for ferromagnetic systems 
based on the Migdal transformation of Ising spins gives 
inaccurate results in high d. However, if one uses the 
parametrization of the 0^ theory, one recovers dc = 4 
within the RSRG Migdal approach even for ferro- 
magnetic systems. 

In any case, notice that our RSRG method does not 
need an explicit geometrical definition of the RG trans- 



formation. Therefore it bypasses the problems related 
to a). In some sense, this is similar to the phenomeno- 
logical RG method |Q where the RG transformation is 
defined implicitly through finite size scaling arguments. 
Remarkably, phenomenological RG calculations are quite 
accurate. 

With respect to point b), as discussed above the ab- 
sence of relevant fields makes truncation errors much less 
important than in ordinary applications of the RG. Fur- 
thermore we note that in the KPZ problem one has to 
compute exponents depending only on the RG transfor- 
mation at the fixed point, i. e. the critical parameter. 
This is profoundly different from what happens in Ising- 
like problems, where some exponents (as, for example, 
the correlation length exponent i^) depend on the deriva- 
tives of the RG transformation around it: As a conse- 
quence i^-type exponents are rather difficult to estimate 
since, even if the location of the fixed point is deter- 
mined accurately, the computation of the derivatives is 
much less precise. No exponents of such type exist in the 
KPZ case. This is, in our opinion, a further reason for 
the great accuracy of the new method with respect to the 
usual RSRG. 

As a final point, it is worth discussing the current 
limitations of the method. It is clear from the results 
presented, that a most important role in the method is 
played by the choice of the parametrization of the scale 
invariant dynamics. This is particularly true since one 
deals with a monoparametric description of the growth 
process: If one could easily introduce several parame- 
ters and study their flow under the RG, the stability of 
the results when details are changed would be greatly 
improved. Within the present framework, the inclusion 
of additional parameters is however not straightforward. 
The problem is that additional RG equations are pro- 
vided only by the use of equation (11) with 6, 6^, and so 
on. This requires the computation of cl;^(6, k) on systems 
whose size becomes quickly prohibitively large. A re- 
markable improvement of the method would therefore be 
the identification of additional RG transformations inde- 
pendent from Eq. (^l]) . Despite this difficulty, we believe 
that the theoretical framework presented here constitutes 
an important new element in the field of surface growth 
and deserves further investigation, in particular with re- 
spect to possible applications to other open problems, 
and generalization to deal with time-dependent proper- 
ties. 

In summary, in this paper we have presented a real 
space renormalization group method developed to deal 
with surface growth processes. The new method over- 
comes the difficulties inherent to standard real space 
renormalization group analysis of anisotropic situations. 
It is based on the definition of anisotropic blocks of 
generic scale and of a parametrized effective dynamics for 
the evolution of such blocks. Imposing the surface width 
to be the same when using different scales (different block 
sizes) we write a renormalization group equation. Its as- 
sociated fixed points define the scale invariant effective 
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dynamics and permit to determine the roughness expo- 
nent a. 

We have employed the new method to study the 
Kardar-Parisi-Zhang and the Edwards- Wilkinson univer- 
sality classes. In particular, for KPZ we compute the a 
exponent in dimensions from d = 1 to rf = 9. The results 
are in very good agreement with the best numerical esti- 
mates in all dimensions. Moreover, we present analytical 
calculation excluding the possibility of KPZ having a fi- 
nite upper critical dimension. On the other hand, well 
known results for the EW universality class are obtained, 
confirming the generality of the method. 
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APPENDIX A: THE RSRG METHOD FOR 
LARGE B 

In this Appendix we investigate the behavior of the 
method for large values of b. At the fixed point x*, 
Eq. (|l|) reads 



uj\b^,x*) = cuj'^{b,x*) + 2uj^{b,x*). 
If we now assume that, for 3> 1, 

oj'^ib, x) ~ 6^" [A{x) + B{x)b-'^ + ... 



as it should, we find that Eq. (Al) becomes 



A{x*)[l - cA{x*)] - 2cA{x*)B{x*)b-"^ 
-2A(a;*)&"^" + subleading terms = 

We see that for 6 — > oo the fixed point tends to 

xlo ■■ cA{xl^) = 1 

whereas for large but finite b 



with 



A = minlct;, 2a} and D{x*) 

B{x*) ifcj<2a, 
A{x*) + B{x*) iftj = 2a, 
A{x*) ifcj>2a. 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 



(A6) 



2\nb 

In + cA{xl) + cB{xl)b-'^ + . . .] 

In [1 + 6-2" - 2cD{x*^)b-^ + cB{x*^)b- 
2hn> 



(A7) 



and converges to the exact value for 6 oo. Note that 
only the finite b corrections depend on x*. In particular, 
if < 2q; 



a = a — 



cB{x^ 
2 In 6 



If w > 2a 



2 In 6 



(A8) 



(A9) 



With respect to the stability of the fixed point one finds 



R'{x* 



(AlO) 



This means that, as it should be expected, the fixed point 
becomes more and more stable as b increases: for larger 
b fewer RG iterations are necessary to reach the scale 
invariant regime. 

These results for large b are not surprising. When 
the systems become large, the effect of the boundary 
conditions is clearly small and also the choice of the 
parametrization of the scale invariant dynamics tends to 
become irrelevant, since parameters that are not included 
in the explicit parametrization are generat ed by th e RG 
procedure on large sytems. The formulae (AS- AlO) cer- 
tify that the RG method is asymptotically correct. 



APPENDIX B: COMPUTATION OF uj^{B,eK). 

In this Appendix we present the derivation of Eq. ([5S|). 
Let us consider b and d arbitrarily large but finite, so that 
Efc < b'^. We have 



6"-! 



1 



71=1 



where 



6''-l 

y - 

^ o 



(Bl) 



Hence 



1 + b<^m} 



Oiel 



(B2) 



(B3) 



The RG estimate d of the exponent a is given by 



and 
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0{el 



(B4) 



The roughness of a system of size b is therefore 



n 



n=l 

yii-y) 



b\k 



dy. 



(B5) 



Fog large b and d and infinitely strong lateral growth 
parameter l/cfc, the set of high sites will form, when n — > 
a d-dimensional hypersphere. Hence ri„ will scale as the 
perimeter of such hypersphere 



Jd-l)/d 



0. 



(B6) 



Similarly for n ^ b'^ the low sites will form a shrinking 
hypersphere and 0„ '-^ {b'^ — ny'^~^^^'^ for n ^ b'^. Hence 
it is reasonable to assume 



ud-1 



(B7) 



with n{y) - for y and n{y) ~ {l-y)(d-^)/d 

for y — > 1. The form of ^(y) for intermediate values of 
y is not known, but we expect it to be nonsingular and 
not dependent on b. In conclusion 



with 



uj\b,ek)=¥'+'ekfi 



Jo n{y) 



(B8) 



(B9) 



a finite geometrical constant. 
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FIG. 1. Covering procedure of a particular microscopic sur- 
face with cells of size Lk x hk. Blocks below (above) the 
surface are considered to be occupied (empty). 
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FIG. 2. Growth rates for the KPZ dynamics in a typical 
configuration. 
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FIG. 4. Results for d = 1 and different Ahmax- The inter- 
sections between solid (representing the 1/2 log2 ^2(2:)) and 
dashed lines (representing the 1/2 log2 [-F4(a;)/F2(a:)]) give x* 
on the horizontal axis and a on the vertical axis (see Eq. (pl|)). 
Observe that, as /S.hmax is increased, the exact value a = 1/2 
is rapidly approached. 
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FIG. 5. Value of a as a function of the dimension obtained 
by the application of the method with small 6, compared with 
numerical data by Ala Nissila el al. The best fit to the RG 
values gives slope —0.7. We expect this exponent to converge 
to —1 for large dimensions. 
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FIG. 7. Plot of the RG estimate for the roughness exponent 
a as a function of the inverse system size. Empty symbols are 
for data obtained with periodic boundary conditions. Full 
symbols refer to open boundaries. 



G- 


- &' b = 


2 


B- 


-Hb = 


4 




-Ob = 


8 




-Ab = 


16 




-^b = 


32 



1000 



0.50 - 



10 



J3- 



a 



^-<1 

-A 

■-[] 



# - -#d = 1 open 
■ - -Hd = 2 open 

♦ - -♦d = 4 open 
G - = 1 closed 
□ -^d = 2 closed 
A - -Ad = 3 closed 
O-Od.A closed 



0.55[i^ ^--^b = 32 - 100 r 



^a:--.;;f-^ • 

-0 f o ®- 



0.45 I ' • — 

0.60 0.70 



I ' ' ' ' ' ' ' ' ' 

0.00 0.10 0.20 0.30 0.40 0.50 

1/b 



" FIG. 8. Value of the fixed point parameter x* as a function 

FIG. 6. Behavior of the curves ai,{x) in d = 1 (see Eq. (^) of the inverse system size. Empty symbols are for data ob- 
and (|37[)). tained with periodic boundary conditions. Full symbols refer 

to open boundaries. 
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FIG. 11. Representation of a surface configuration with 8 
"high" sites and fc"* — 8 "low" sites. 
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FIG. 9. Value of the exponent a computed using the 
parametrization Eq. ( po| ) for the KPZ dynamics for various 
values of 7 in d = 2 and d — 4. 
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FIG. 12. Comparison of the form of the function Qn com- 
puted with a:: = 0, used in the analytical calculation (solid 
line), with the function determined numerically in d = 7 for 
X = oo (dashed line). For intermediate values of x, the func- 
tion fin lies between the two curves plotted. 
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FIG. 10. Value of the exponent a computed using the 
parametrization Eq. ( ]4l| ) for KPZ dynamics for various values 
of 7 in d = 2 and d = 4. 
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d 

FIG. 14. Comparison of the dependence of a on d com- 
puted with the RG in the large-d hmit Eq. ( [53[ ) with the 
numerical data of Ala-Nissila et at. In the inset, plot of 
the value of the fixed point parameter x* computed numeri- 
cally for small d and compared with the analytical expression 
Eq. (H. 



d 


1 2 3 


4 


5 


6 


7 


8 


9 


X* 
Q!RG 
dnum 


0.726 2.77 6.96 
0.507 0.363 0.294 
0.5 0.387 0.305 


15.91 
0.238 
0.261 


31.96 
0.206 
0.193 


63.5 
0.182 
0.18 


124.5 
0.164 
0.15 


242 
0.149 


468 
0.137 



TABLE I. Values of the fixed point parameter x* and of 
the roughness exponent computed with the RG with 6 = 2, 
compared with the numerical results of Ala-Nissila et al. 
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